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Abstract 



The gravitational radiation from compact pulsar-like stars depends on the state of 
dense matter at supranuclear densities, i.e., the nature of pulsar (e.g., either normal 
neutron stars or quark stars). The solid quark star model is focused for the nature 
of pulsar-like compact objects. Possible gravitational emission of quark stars (either 
fluid or solid) during the birth and later lifetime is discussed. Several observational 
features to distinguish various models for pulsar-like stars are proposed. It is sug- 
gested that the gravitational wave behaviors should be mass-dependent. Based on 
the data from the second LIGO science run, the upper limits of R ■ 9^^^ and thus 
mass, R: radius, and 6: wobble angle) are provided for those targets 
of millisecond pulsars. 
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1 Introduction 



On the one hand, it is now a very critical time for detecting directly grav- 
itational waves. Despite the negative detection results, some limits on the 
physical parameters of pulsar-like star's interiors can be obtained with the- 
oretical calculations which relate the waves to star's nature in literatures, 
base d on the experiment data of gravitational waves. Andersson & Kokko- 
tas ( Andersson fc Kokkotasl . 199^ calculated the eigenfrequencies of the /, 



p, and w modes, and four id that the frequencies depend on star's mass and 



radius. Numerical results (jBenhar. Berti fc FerrarL „199ft) showed that the in- 



formation of neutron star structure could be carried by gravitational waves of 
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the axial tf-modes. The inverse problem for pulsating neutron stars, based on 
the studies of thos e modes, was discussed extensively by Kok kotas, Aposto- 
latos & Andersson (jKokkotas. Apostolatos &: Anderssonl . 120011 1 . The frequen- 
cies and damping behaviors were r eexamined for different equ ations of state 
(iBenhar. Ferrari &: Gualtieri . 2004 ), and recently, Tsui & Leung ( Tsui &: Leung 
200,^^ applied an inversion scheme to determine star's mass, radius, and den- 
sity distribution via gravitational wave asteroseismology. From the theoretical 
work above, one sees then that a positive result of recording gravitational wave 
signal should improve significantly the knowledge of matter at supranuclear 
density ^ . 



On the other hand, pulsar-like compact stars were discovered since 1968, but, 
unfortunately, we are still not sure about their nature. It is conventionally 
believed that these compact stars are normal neutron stars. However one 
can not rule out the possibility, that they are actually quark stars com- 
posed of quark matter with de-confined quarks and glu ons, either from the 
first principles or according to their vario us observations ( Glendennine . 2000l : 
Lattimer fc Prakash . 2004; Weber . 2005[ ). Quarks are fundamental Fermions 
in the standard model of particle physics, and the existence of quark mat- 
ter is a direct consequence of the asymptotically free nature of the strong 
interaction, which was proved in the regime of non-Abelian gauge theories 
and confirmed by high-energy collider experiments. In this sense, to affirm 
or negate the existence of quark stars is of significantly fundamental mean- 
ing for understanding the nature's elementary strong interaction. Quark stars 
with strangeness are popularly discussed in literatures, which are called as 
strange (quark) stars, whereas nonstrange quark stars would also be possible 
if nonstrange quark matter could also b e stable at zero pressure due to color 
confinement (|Menezes fc Melrosd . 2005| ). We do not differentiate between the 
terms of quark stars and strange stars here, supposing both kinds of the stars 
have a quark surface to confine quarks and gluons. 



Therefore, one of the key points for today's astrophysicists is understanding 
the nature of pulsars: to find competitive evidence for quark stars or for normal 
neutron stars. Effective methods to do depend on (i) the minimu m spin per i- 
ods, (ii) the mass-radius relations, and (iii) the surface differences (2£u, 2003bl ). 
However, since they are relativistic, pulsar-like compact stars could be poten- 
tially gravitational wave radiators. The features and strength of this emission 
would reflect their nature, and one may probe into pulsar's interior with grav- 
itational waves. Recently, the Laser Interferometri c Gravitational wav e Ob- 



servatory (LIGO) has monitored 28 radio pulsars ( Abbott et al. . 2005[ ). and 



obtained upper limits of gravitational radiative intensities for the pulsars. If 



^ It may be worth noting, however, that gravi tational wave or graviton might not 
exist at all in the de Sitter Universe ( LiuL 2004h if gravitational wave would still not 
be detected in the next decades. 
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pulsars are quark stars, how and what can one constraint the pulsar's nature 
based on the data? We are trying to analyze relevant issues in this paper. 

There are actually two motivations to study gravitational waves using pulsars: 
(i) as a tool and (ii) as a source. These two result in different physics. The 
radio pulsar timing array depends on the "clock" nature, not on the internal 
structure of pulsars, while the direct detection of the waves from isolated 
and/or binary pulsars depends on the internal structure of pulsars. The former 
is not discussed in this paper. Certainly, a successful detection of gravitational 
waves from pulsar-like stars is very important for understanding (i) the physics 
of matter at supranuclear density and (ii) the physics of gravity. 



2 Quark stars during their birth 



A mass formula for strangelet analogous to the Bethe -Weizsacher semi-empirica l 
mass function in nuclear physics was introduced bv lDesai. Crawford fc Shaw 
Neglecting the symmetry and Coulomb energy terms, which are not 
important in the case of strange quark stars, one has E = + Eg, where the 

volume and surface energy contributions are E^ = e^A' and E^ ~ O .lenA''^^^, re- 

spect ively, with eo = (880 ~ 890) MeV and A' the baryon number ()Desai. Crawford fc Shaw . 
1993| ). By introducing baryon number density rib = A'/(AitR^/3) {R the ra- 



dius of spherical quark matter) and the surface energy per unit area 7, we 
have R = (3A'/(4:7m h)V^^ and thus 7 ~ 0.02eon/^ ~ 20 MeV/fm^. However, 
Bagchi et al. obtained a surface energy (or surface tension), 7, which 



varies from about 10 to 140 MeV fm~^, depending on stellar radius. It is worth 
noting that the surface energy, which will be neglected in this paper, due to 
an electrostatic field of <~ 10^^ V/cm on a star's surface is several orders 
smaller that 10 MeV fm~^. We consider thus this poorly known parameter 7 
to range between ~ 10^ MeV fm'^ and ~ 10^ MeV fm'^. 



2.1 Quark nuggets ejected 



A hot turbulent bare quark star may eject (or "evaporate") low-mass quark 
matter (quark nuggets). Phenomenologically, there could be 4 steps to create 
a quark nugget (Fig.l). Two parameters, which are not know with certainty 
yet, describe the turbulent nature of protostrange stars: the conv ective scale L 
and velocity Vt- It is possible that 4 ~ lO^'"^ cm, Vt ~ 10^ cm/s ( Xu fc Bussd . 



200l[ ). In step (a), the typical kinematic energy of turbulent part is Et ~ pl^v^ 



while the surface and the gravitational energies, which may prevent the fluid 
part to eject, are Eg ~ 7/^ and Eg ~ GmMr / R"^ (quark star's mass M and 
its radius R, lump's mass m and radius r), respectively. The height of a bulge 
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(a) (b) (c) (d) 

Fig. 1. Forth steps to eject a quark nugget during the formation of a quark star. 

increases if + Eg < E^, as is shown in step (b). The contact area between 
the star and the bump becomes smaller and smaller, as shown in step (c), 
because of (i) a limited scale of kinematic flow with order, and (ii) a total 
momentum which tends to separate nugget from star. A quark nugget form 
finally, with a velocity of V^n <~ ^^t, step (d). The density of a bare st range 



star with mass <~ Mq is nearly uniform (|Alcock. Far hi Sz 01intol . ll98fih . and 
the mass of such a star can be well approximated by. 



M ~ -rf(45). 



where the bag constant 5 = (60 ~ 110) MeV/fm^ i.e., (1.07 ~ 1.96) x 10^^ 
g/cm^. For simplicity, we apply this approximation throughout this paper. 

What kind of quark nugget could be ejected? Let's compare the energies below, 

Et R'^v? 1 tV? Ef Pnv?r s i •> 



where po — 4 x 10^ g/cm^ is the average density of quark stars with mass 
<~ Mq, R = R6X 10^ cm, 7 = 7100 x 100 MeV/fm^. Both Et/Eg and Et/E^ 
should be much larger than 1 if quark nuggets are ejected. We conclude then 
from Eq.(2) that: (i) it needs Vt > 10^ cm/s that nuggets with mass > 10^° 
g can be ejected; (ii) nuggets with higher mass could be ejected from nascent 
quark stars with lower mass; (iii) ejection of nuggets with planet masses (e.g., 
the Earth's mass ~ 10^^ g) could be possible if Vt > 10^ cm/s; (iv) a nascent 
quark star with weak turbulent (e.g., Vt < 10^ cm/s) could hardly evaporate 
nuggets with low mass (e.g., < 10^^ g). Hadrons may not evaporate from the 
surface of a quark star with temperature being much lower than 100 MeV if 7 
is order of 100 MeV/fm^. The unknown parameter Vt would be probably much 
larger than 10^ cm/s since observations show that pulsar-like stars receive a 
large kick velocity (of order a few hundred to a thousand km/s) at birth. 
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These ejecta could be captured by the center strange star (with mass M and 
radius R) for Vt should be smaller than the escape velocity J2GM/R ~ 10^° 



cm/s. We then conclude that quark planets (i.e., planets of quark matter, 
with possible mass from much lower than earth mass to about a Jupiter mass) 
would form simultaneously during the formation of a strange star with strong 
turbulence if the surface energy is reasonable (7 < 10^ MeV fm~'^). This 
consequence could be tested by searching quark planets around pulsars. This 
is also an alternative mechanisr n for creating pulsar-pl anet systems, the first 



5ar-pl; 

one of which was discovered by IWolszczan fc Frail |l992lf ) . 



2.2 Global oscillation of quark stars 



A spherical fluid should be deform by rotation. For a quark star of incompress- 
ible fluid with axis-symmetric deformation, the star's surface can be described 

as, 

r(e) = Ai?[l-X:a,P,(cos^)], (3) 

i=2 



where the parameter A is determined by a constant volume of AnR^/S. In case 
of 0:2 ^ 1 (and thus terms with orders being higher than are neglected), 
the quadrupole deformation {i = 2), which will only be noted as following, 
is actually of an ellipsoidal figure with an ellipticity e = a2 and axes a = 
R/y/T^ = b> c = R{l-e). 

In case of a rotating fluid star with T/|iy| <^ 0.1375, the ellipticity e = 
(/ — /o)//o and the eccentricity e = 1 — c^/a^ are related by e ~ e^/3, where 
/(fi), as a function of the spin frequency Q, is the total moment of inertia, with 
Jq = 2Mi?^/5 for stars with uniform density. The total energy, E, of a rotating 
fluid star is then the sum of gravitation, rotation, volume, and surface energies . 



A cal culation similar to the nuclear liquid drop model ^ (jGreiner fc Maruhn 
h99(h shows 



E = -Egravi + + Eg + E^ot 

= £„+(A, + ^.).^ + ^^. (4) 

where Eq is the energy for a non-rotating spherical star, L = IQ the star's 
angular momentum, I{e) the momentum of inertia, and the coefficients Ag 

^ The Coulomb interaction there is replaced by gravitational interaction, with a 
sign-change. 
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and As measure the increases of gravitation and surface energies, respectively, 
of the star. 



The surface energy contribution to the deformation should be negligible since 
the parameter 



is much smaller than 1 unless the star is very low-massive (i.e., mass <~ 10^^ 
g for strangelets) and/or 7 ^ 100 MeV/fm^. 

By using d-E'(e)/d£: = with Eq.(4), neglecting As, one comes to the rotating 
figures of conventional Maclaurin spheroids ^ , 

n ^ or, eo{P) ^ ^ ^ 3 x lO-^PfoL, (7) 



where cq and Eq are the eccentricity and the ellipticity, respectively, of a fluid 
star at equilibrium, which are only a function of spin, provided that po is a 
constant, the spin period P = 27t/Q = Pioms x 10 ms. One can also have 

2(A, + A,) + — ^>0, (8) 



d£2 V 3 ' ' /o(l + £) 



which means the equilibrium state of e = eo is stable. We expect that a fluid 
quark star, either at birth or in its lifetime, should oscillate globally around Eq, 
though the oscillation frequency and amplitude are still not computed now. 
If the amplitude is large, this quadrupole oscillation could result in enough 
gravitational radiation to be detected by future advanced facilities. 



Additionally, this oscillation modes may help us to distinguish fluid and solid 
quark star models. A protoquark st ar should be in a f luid state, but the star 

Therefore, we should 



may solidified soon after its birth (jXuL l2003a . 
find ways to differentiate solid and fluid quark stars in their later time. Due 
to the strong shear force, a solid quark star may oscillate with much smaller 
amplitude but much higher frequency than that of a fluid quark star. Besides 
detecting the gravitational wave, pulsar timing in radio, optical, as well as 



Note that the approximation of eo{P) can only be valid for Pims 3> 1- 
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X-ray bands could be effective if the star has pulsed emission, because the 
rotation period changes during the oscillation. 



Could the oscillation amplitude be larger enough to fragmentate a fluid quark 
star? Unlike the case of nuclear fission, where the Coulomb interaction favors, 
the gravitational interaction prevent a quark star to fragmentate. It needs 
then energy to "exc ite" the star, with order of Eg^ ~ 0.6G(M/2)^/i?. If pulsar 
kick (|Popovl . 12004^ . with energy of i?ki ~ 0.5MV^, plays a major role to 
fission a quark star, we note that the kick velocity Vk should be greater than 



~ -Ry27rG'po/5 ~ 6 x lO^R cm/s. Therefore, a quark star with large kick (or 
turbulent) energy and low mass could fragmentate. 



3 Gravitational waves due to r-mode instability 



In Newtonian theory, a rapidly rotating fluid Maclaurin spheroid is secularly 
unstable to become a Jacobi spheroid, which is non-axisymmetric, if the ratio 
of the rotational kinetic energy to the absolute value of the gravitational poten- 
tial energy T/IVTI > -137 5. In the general relativisti c case, as was shown by 
Chandrasekhar ( 197Cll l and Friedman fc Schut3 ra.gravitafoual radiation 
reaction amplifies an oscillation mode, and it is then found that the critical 
value of T/|W^| for the onset of the instability could be much smaller than 
0.1375 for neutron stars with mass >~ Mq. This sort of non-axisymmetric 
stellar oscillations will inevitably result in gravitational wave radiation. 



A kind of oscilla tion mode, so call e d r-mode, is focused on in the li teratures 

since t he work by Andersson ( 1998| ) , Friedman fc Morsink ( 1998| ) , and Lindblom. Owen fc Morsink 
f 1998[ l. The r-mode oscillation is also called as the Rossby waves that are ob- 
served in the Earth's ocean and atmosphere, the restoring force of which is the 
Coriolis force. This instability may increase forever if no dissipation occurs. 
Therefore, whether the instability can appear and how much the oscillation 
amplitude is depend on the interior structure of pulsars, which is a tremen- 
dously complicated issue in supranuclear physics. 



To make sense of the generic nature for r-mode instability in different star- 
modes but avoiding an uncertainty in microphysics, we just present a rough 
calculation below. The critical angular frequency, Q, limited by the gravita- 
tional radiation due to r-mode instability, are determined by comparing the 
damping/growth times due to the various mechanisms for energy dissipation 
in a given star-model. 



1 ^ 1 _^ 1 

Tgw T~^v Toy 



0, 



(9) 
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where the timescales for the instabihty are estimated to be ( Madsen . 1998| ). 



'gw 
_ns 



sv 

_ss 
'bv 



-3.85 X lO^^Q'^ M-^R-^, 
1.01 X IQ-^M-^R^T^, 



t{^^ = 1.29 X lO^^Q"^ M-^R^T-^, 



(10) 



= 1.85 X 10-9«.5/3M-5/9i?ii/32-5/3^ 



5.75 X IQ-^m^oln^R^T-^, 



where superscript "ns" ("ss") denotes neutron (bare strange) star model, Tsv 
and Tbv are the dissipation timescales due to shear and bulk viscosities, respec- 
tively, and as the coupling constant of strong interaction, T the temperature, 
mioo the strange quark mass in 100 MeV. It is assumed in Eq.(lO) that the 
modified URCA process dominate the weak interaction in normal neutron 
stars. 

Strange stars could have very low masses, even of 1O~^M0 (a strange star 
with planet mass could also be called as a strange guarA; planet). A low-mass 
normal neutron star could not be lik ely because of (i) the minim um mass 
of a stable neutron star is ~ O.IM© ( Shapiro fc Teukolskvl . Il983| ). and (ii) 
the gravitation-released energy of possible low-mass neutron star during an 
iron-core collapse supernova. Eg ~ GM'^/R, should be much smaller than 
~ ]^q53 gj,g (-^^g ^Yie relation M ~ R~^ for low-mass neutron stars, but it 
seems no t successful i n modern supernova simulations even in case of Eg ~ 
Iq53 gj,g ( Buras et al. . 20031 : Liebendoerfei . 2004). Therefore, in the following 
calculations, we consider low mass issues for bare strange stars, but apply only 
canonical mass of 1.4Mq for normal neutron stars. 

No r-model instability occurs in a solid star since the Coriolis-restoration force 
does not work here. The break frequency of low-mass bare strange stars could 
be approximately a constant. 



IGM 



1.1 X 10^ s" 



:ii) 



with a pref actor of ~ 0.65 at most for M ~ Mq and i? ~ 10^ cm ( Glendennind . 
2000f ). where the bag constant B = QO MeV/fm^. Let's express the critical 



frequency in unit of Qq, through calculations based on Eq.(9). It is found in 
Fig. 2 that: (i) gravitational wave radiates more likely from proto-neutron stars 
than from proto-strange stars; (ii) the r-mode instability could not occur in 
fluid bare strange stars with radii being smaller than ~ 5 km (or mass of a 
few O.IM0) unless these stars rotates faster than the break frequency; and (iii) 
the above conclusions do not change significantly in the reasonable parameter- 
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6 7 8 9 10 
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Fig. 2. Temperature dependence of the critical angular frequency, in unit of Qq, due 
to r-mode instability for normal neutron stars and bare strange stars. The critical 
lines are numbered from "0" to "6". "0": normal neutron stars; "1": bare strange 
stars (BSSs) with mass M = Mq, the bag constant B = 60 MeV/fm~^, the strange 
quark mass rus = 100 MeV, and the coupling constant Og = 0.1; "2": BSSs with 
M = O.IM0, B = 60 MeV/fm-3, = 100 MeV, and = 0.01; "3": BSSs with 
M = O.IMq, B = 60 MeV/fm-3, = 100 MeV, and as = 0.1; "4": BSSs with 
M = O.OOIM0, B = 60 MeV/fm-3, ms = 100 MeV, and as = 0.1; "5": BSSs with 
M = Mq, B = 110 MeV/fm-3, rris = 100 MeV, and as = 0.1; "6": BSSs with 
M = Mq, B = 60 MeV/fm-3, nis = 100 MeV, and = 0.9. 



space of B, as, and mg. To estimate strength of gravitational waves from this 
oscillation, one needs to simulate the nonlinear increase of the instability, 
which is not certain yet. 



Some recent observation s in X-ray a stronomy could hint the existence of low- 
mass bare strange stars ()XuL l2005a ). The radiation radii (of, e.g., IE 1207.4- 
5209 and RX J1856. 5-3754) are only a few kilometers. No gravitational wave 
emission could be detected from such fluid stars even they spin only with a 
period of ~ 1 ms. 
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4 Gravitational waves from wobbling stars 



Gravitational wave radiation from pulsars could be classified as (i) emission 
due to the normal modes of oscillation of fluid matter (e.g., the r— mode 
instability discussed in §3), and (ii) emission due to solid deforming which 
will be focused in this section. Protoquark stars should be in a fluid state 
when their temperatures are order oflOMeV, but would be solidified as they 



cool to very low temperatures (j2£u|, |2Q03a|) . Assuming the ellipticity of a solid 



quark star, with an initial spin period Pq, keeps the same as that of the star 
just in its fluid phase, we expect realistic ellipticity e of the star in the solid 
state satisfies 

£o(Po) > e > eo{P) (12) 



as the star spins down to a period of P, due to the shear force that prevent the 
star to deform. For the Earth which is spinning down, the mean density is 5.5 
g/cm^, the ellipticity observed e = 0.00223, one has eo{2A hours)= 0.00287 >~ 
e from Eq.(7). Additionally, a free precession mode called the Chandler wobble, 
which was formulated and expected by Leonhard Euler, with a period of 435 
days suggests a good approximation of e too. This means that the suggestion 
of Eq.(12) is at least workable for the equilibrium figure of the Earth ^ . 

However, strain energy develops when a solid quark star spins down, which is 
proportional to [eo{Po) — e]^. The first glitch of the quark star occurs when the 



strain energy reaches a critical value (jZhou et al.l . l2004 ). Subsequent glitches 



take place too as the stellar stress increases to critical points. We note that the 
star's ellipticity decreases after glitches due to (i) the release of strain energy 
and (ii) possible increase of mean density if stellar volume shrinks as the star 
cools. 

A star's ellipticity could be observed if the star's free precession mode is discov- 
ered due to the relation e ~ P/ Pprece, where Pprece is the pre cession period. Evi- 



dence for free precession for PSR B1828-11 was provided bv lStairs. Lvne fc Shemar 



Deriodic timing variation and changes in pulse duration. The derived 



ellipticity is ~ P/Pprece — 0.405 s/500 d ~ 10~® if the precession is really in 
a free mode. But, based on Eq.(7), the ellipticity is eo(405 ms) ~ 2 x 10^^, 
which is 2 orders larger than that expected. This conclusion conflicts with 
Eq.(12). Why is PSR B1828-11 so "round"? An answer proposed here is that 
this observed precession mode could actually be torqued by external force, 
rather than free. In fact, a disk-forced precession model for PSR B1828-11 



^ Among solid astrophysical objects, only the Earth is investigated mostly exten- 
sively. One has then to learn from solid geophysicists in order to develop a compre- 
hensive astrophysical model for solid quark stars. 
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was provided bv lOiao et all (|2nn3h . but other source torques by, e.g., magne- 
todipole radiation or planet (s), are also possible. The precession period could 
be very long (and thus results in a much small ellipticity if assuming free 
precession) if torques are not strong. Due to the torque exerted by the Moon 
and the Sun, the spinning Earth shows a forced precession period of 25,700 
yr, being much larger than 435 d, the period of free precession. By the way, 
planet-torque induced precession modes could be effective to search planets 
around pulsars. Nevertheless, a free precession mode could have been dis- 
covered in a bursting r adio source, GCR T J1745-3009, in the direction of the 
Galactic center region ( Zhu Sz Xu . 2005[ ). The eUipticity derived from P/Pp 



prece 

is ~ 10/77 ~ 2 X 10~^Pioms- Based on Eqs.(7) and (12), the spin period should 
be P > 0.1 s, which is very reasonable for normal pulsars. If the pulsar's el- 
lipticity tracks almost eo{P) due to glitches, e >~ So{P), then its ellipticity 
e >~ 2 X 10~^ and period P >~ 0.1 s. Statistically, a pulsar with P >~ 0.1 
s could be a radio nulling (even extremely, with null fraction > 9 0%) pulsar 
( Biggj . ll992h . which fits the model provided by (|Zhu fc XilboO,^ . 



A pulsar must be non-axisymmetric in order to radiate gravitationally. A 
wobbling pulsar, either freely or forcedly, may thus radiation gravitational 
waves. This wave results in a perturbed metric /t^„^ of space-time (q,,,^ = p,,,,, + 



hii.„), which is order of /tn being given by ( Jaranowski. Krolak fc Schut d . 1998^ 
Jones fc Andersson . 2002f l. 



ho 



15c4 



dP^ 



^ ^ 2.8 X 10 



t-f^6"kpc-' lOms^^' 



(13) 



where approximations / ~ 0.4Mi?^ and M ~ AnR^po/S are applied for solid 
quark stars in the right equation, the pulsar's distance to earth is ci = (ikpc x 1 
kpc, 6 is the wobble angle. LIGO is sensitive to hight frequency waves, which 
recently puts upp er limits on ho for 2 8 known pulsars through the second 
LIGO science run ( Abbott et al. . l2005[ ). The upper limits are order of 10~^^, 
which means approximately an limit of ei?g < 10~^. This observation may not 
conflict with Eqs.(7) and (12) if the pulsars' radii are slightly smaller than 
10 km. In fact, poss i ble can didates for low- mass quark stars with small radii 



ar e proposed bylXul (l2005a[ l. Only three normal pulsars are listed in Table 1 
of lAbbott et"d1 ()200,5t ): the Crab (B0531+21), J1913+1011, and B1951+32, 



others are millisecond pulsars. 

It is conventionally suggested that millisecond pulsars are recycled, with a 
spinup history during binary accretion, and the ellipticity may not follow 
Eq.(12) since eo{P) > e for such stars during their spinup phases. Nonethe- 
less, since millisecond pulsars are now in a spindown phase, we may simply 
suggest eo{P) ~ e, which is valid for the Earth, and calculate the upper 
limits of pulsars' radii and thus masses according to Eqs.(7) and (13). Note 
that Eq.(12) is still effective if millisecond pulsars are born during an AIC 
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(accretion-induced collapse of white dwarfs) process (j2£u|, |2Q05a|) . The upper 
limits of pulsars' radii and masses are calculated, which are presented in Ta- 
ble 1. Large radius limits (> 10 km; and thus high mass limit, > Mq) are for 
the three normal pulsars, which does not show any constrain on the equation 
of state since pulsar's maximum mass is only ~ 2Mq. However, the LIGO's 
observation suggests a small limit of R ■ 9^^^ (a few kilometers), and thus a 
small limit of M ■ 9^^^ (between 10~'^Mq and 10~^ Mq), for those millisecond 
pulsars. This hints either a low mass (<C ^o) of pulsars or a small wobble 
angle 9. Factually, in order to explain its polarization behavior of radio pulses 
and the integrated profile (pulse widths of main-pulse and inter-pulse, and 
the separation between them), the fastest^ rotating millisecond pulsar PSR 
J19 39-F2134 (B1937 +21) is supposed to have mass < 0.2Mq and radius < 1 
km (IXu et alilimh . This conclusion does not conflict with the value in Table 



1. Observationally, the precession angle (i.e., the angle between the spin vector 
and total angular momentum), which is smaller than 9, should be smaller than 
the pulsar radio beam angle (a few tens of degrees). The radius upper limits for 
the millisecond pulsars could be about two times the values listed in Table 1, 
and then the mass limit of (3 ~ 10) times, if 6' ~ (1° ~ 10°). Actually, a small 
radius of pulsar-like stars is not surprising since we have detected six Cen- 
tral Compact Obje cts in supernova remnant s (CCO s) with black-body radius 
of (0.3 ~ 2.4) km (IPavlov. Sanwal fc Teteil . EoJ, and seven Dim Thermal 
Neutron stars (DTNs) with possible radius of a few kilometers f Haberl 2005| ). 



Elliptic deformation v.s. bumpy distortion. The lack of symmetry about a 
pulsar's rotation axis could be the result of either elliptic deformation because 
of rotation, which is discussed in §2 and §4. or bumpy distortion (i.e., localized 
mountains on stellar surface). Strong magnetic field is suggested for producing 
the bumps, since the magnetic and rotating axes are generally not aligned. 
However for quark stars, this mechanisms may not work due to a relatively 
negligible magnetic force, {B"^ / 8n) / (pqc^) ~ 10~^i?i6, even for fields {B = 
BiQ X 10^^ G) as strong as 10^^ G. G litches of sol id quark stars could produce 
bumps, with a maximum ellipticity ()Owen . 2005[ ). 



lO-H^^)RAl + 0.08ARlr\ (14) 



where amax is the stellar break strain. This ellipticity is larger for low-mass 
quark stars due to weaker gravity. But for normal neutron stars, Owen ( 2005h 
derived the maximum elastic deformation, emax, induced from shear stresses 
is typically only 6.0 x 10~^. The real ellipticity of a quark star could be far 
from emax, but may be app roximatelv eri(P), due to stress releases through 



star-quake induced glitches (jZhou et al.l . 120041 ) . 



" A 1.39 ms radio pulsar could have been discovered in the globular cluster Terzan 
5 by the Robert C. Green Bank Telescope (GBT). 
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pulsar 


period 


distance 


ho 


radius-y ' 


mass-t^^ 


name 


P (ms) 


d (kpc) 


(10-24) 


(km) 


(Mq) 


B0531+21* 


33.08 


2.00 


41 


19.6 


6.4 


J1913+1011* 


35.91 


4.48 


51 


18.6 


5.5 


B1951+32* 


39.53 


2.50 


48 


22.4 


9.4 


B0021-72C 


5.76 


4.80 


4.3 


2.6 


0.015 


B0021-72D 


5.36 


4.80 


4.1 


2.4 


0.012 


B0021-72F 


2.62 


4.80 


7.2 


1.5 


0.0030 


B0021-72G 


4.04 


4.80 


4.1 


1.9 


0.0061 


B0021-72L 


4.35 


4.80 


2.9 


1.9 


0.0059 


B0021-72M 


3.68 


4.80 


3.3 


1.7 


0.0043 


B0021-72N 


3.05 


4.80 


4.0 


1.5 


0.0031 


J0030+0451 


4.87 


0.23 


3.8 


4.1 


0.056 


J0711-6830 


5.49 


1.04 


2.4 


3.0 


0.023 


J1024-0719 


5.16 


0.35 


3.9 


3.9 


0.051 


B1516+02A 


5.55 


7.80 


3.6 


2.2 


0.0090 


J1629-6902 


6.00 


1.36 


2.3 


3.0 


0.024 


J1721-2457 


3.50 


1.56 


4.0 


2.1 


0.0084 


J1730-2304 


8.12 


0.51 


3.1 


5.0 


0.11 


J1744-1134 


4.07 


0.36 


5.9 


3.5 


0.037 


J1748-2446C 


8.44 


8.70 


3.1 


2.9 


0.021 


B1820-30A 


5.44 


7.90 


4.2 


2.2 


0.0094 


B1821-24 


3.05 


4.90 


5.6 


1.6 


0.0037 


J1910-5959B 


8.36 


4.00 


2.4 


3.2 


0.028 


J1910-5959C 


5.28 


4.00 


3.3 


2.4 


0.011 


J1910-5959D 


9.04 


4.00 


1.7 


3.2 


0.028 


J1910-5959E 


4.57 


4.00 


7.5 


2.5 


0.013 


J1939+2134 


1.56 


3.60 


13 


1.2 


0.0015 


J2124-3358 


4.93 


0.25 


3.1 


3.9 


0.048 


J2322+2057 


4.81 


0.78 


4.1 


3.2 


0.027 


Table 1 



The upper limits of the radii and masses of 28 pulsars targeted in the sec- 
ond science run of LIGO, which depend on the wobble angle 9. Normal 
pulsars liste d are starred (*), others are millisecond pulsars. The /iq lim- 
its are from lAbbott et al.l Eol), while the pulsar data are obtained through 
ihttp:/ /www.atnf.csiro.au/research/pu ls^gf ps rcat / \ 



5 Conclusions 



Detection of gravitational waves can certainly test the general theory of rela- 
tivity, and also open a new window for us to observe astrophysical phenomena. 
Observations in gravitational wave band may reveal the nature of pulsar-like 
stars, and help to answer numerous questions: Are pulsar-like stars normal 
neutron stars or quark stars? Is quark matter with hight density but low tem- 
perature in a solid state? Do quark stars with low masses exist in the Universe? 
The paper is summarized as follows. 

(1) . During the birth of quark stars, quark nuggets may be ejected, and the 
global quadrupole oscillation around elhpticity eo{P) [see Eq.(7)] could also 
be a source for gravitational waves. 

(2) . No gravitational wave originated from r— mode instability could be de- 
tected for pulsars to be quark stars with low masses or in a solid state, for 
isolated as well as accreting compact stars, even the compact stars spin at 
sub-millisecond periods. 

(3) . Through the second LIGO science run, we may find upper limits of masses 
and radii for millisecond pulsars, although one can not constrain on the masses 
and radii of normal neutron stars. The radius of PSR B1937+21 could be 
smaller than ~ 2 km if its wobble angle 6 is between 1° and 10°. Future LIGO 
observations could put tight constraint on and thus if pulsar-like 
stars are actually solid quark stars. 

(4) . We suggest the precession mode discovered in PSR B1828-11 is not free, 
but could be forced by fossil disk or planets. The geometrical parameters of 
torqued precession might be obtained by precise timing of radio pulses. The 
bursting radio source, GCRT J1745-3009, could be a freely precessing radio 
pulsar with spin period ~ 0.1 s, which could be tested by searching pulsed 
emission from the source. 

As the mass-radius relation for low-mass quark stars (M oc R^) is in strik- 
ing contrast to those of normal neutron stars with low masses (M oc R^^), 
gravitational radiation from low-mass neutron stars and quark stars should 
be very different for processes of both fluid instability and stellar wobbling. 
We have just studied mass-dependent waves for quark stars, but it is also 
necessary to obtain the mass-dependent behavior for normal neutron stars, 
in order to put theoretical models directly in front of gravitational wave ob- 
servations. Additionally, similar to the case of double white dwarf binaries 
()Stroeer. Vecchio fc Nelemansl . boosh . detecting gravitational wave via LISA 
(the Laser Interferometer Space Antenna) from binary quark stars is also 
valuable to constrain the masses and radius. The radiation reaction, tides, 
and mass transfer are very different for binary low-mass quark stars, white 
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dwarfs, and normal neutron stars. 
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